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Abstract 

Extending our earlier treatments of vr°,r/ c and rjb, we study the rj—rj' sys- 
tem and its 77 decays using a model which is a leading version of the con- 
sistently coupled Schwinger-Dyson (SD) and Bethe-Salpeter (BS) approach. 
The electromagnetic interactions are incorporated through a (generalized) im- 
pulse approximation consistent with this bound-state approach, so that the 
Ward-Takahashi identities of QED are preserved when quarks are dynamically 
dressed. To overcome some of the limitations due to the ladder approximation, 
we introduce a minimal extension to the bound-state approach employed, so 
that the Ua(1) problem is avoided. Pointing out which of our predictions hold 
in the coupled SD-BS approach in general, and which are the consequences of 
the specific, chosen model, we present the results for the axial-current decay 
constants of rjs, rjo, and of their physical combinations 77 and 7/, the results 
for the 77-decay constants of 770 and rjs, for the two-photon decay widths of 
r\ and 7/, and for the mixing-independent i?-ratio constructed from them. 
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1. INTRODUCTION 



A particularly interesting example of the applications of Schwinger-Dyson equations to 
hadronic physics (reviewed in, e.g., 0,3), is the approach through consistently coupled 
Schwinger-Dyson (SD) equations for quark propagators and Bethe-Salpeter (BS) equations 
for bound states of quarks. Among various studies of this kind, those of Jain and Munczek 
are judged by many as "the most extensive and phenomenologically successful spectro- 
scopic studies in the rainbow-ladder approximation" || and therefore often chosen |1],|2|,|^-|8| 
as a representative, paradigmatic example of such studies. The essence of such a treatment 
of qq bound states is the solving of the ladder Schwinger-Dyson (SD) equation for the dressed 
quark propagator S(q), and then solving in the consistent approximation, with this result- 
ing dressed quark propagator and with the same interaction kernel, the Bethe-Salpeter (BS) 
relativistic bound-state equation for a qq meson. This procedure is crucial for obtaining 
the mesons from the light pseudoscalar octet as Goldstone bosons when the chiral symme- 
try is spontaneously broken. Thanks to this, a coupled SD-BS approach (notably, || [| ) 
can reproduce the correct chiral limit behavior (crucial in the light sector) simultaneously 
with the realistic results for heavy mesons. In [Q-|5[], the interaction kernel is given by a 
modeled gluon propagator consisting of: a) the well-known perturbative part, reproducing 
correctly the ultraviolet (UV) asymptotic behavior unambiguously required by QCD in its 
high-energy, perturbative regime, and b) the nonperturbative part, which should describe 
the infrared (IR) behavior. Since the IR behavior of QCD is still more or less unknown, this 
latter, nonperturbative part of the gluon propagator is modeled. In |3|-|5| , several forms for 
this IR-part have been used and their parameters varied, with the outcome that results are 
not very sensitive to such variations. Jain and Munczek |3|-|5| have succeeded in reproducing 
the leptonic decay constants of pseudoscalar mesons, and, even more importantly, a very 
large part of the meson spectrum, except for such elusive cases as the 77-7/ system. 

Such an up till now successful and reputable referent model should be tested further by 
calculating other quantities (e.g., electromagnetic processes) to see how well it will do. This 
was our motivation for calculating ir°,T] c , rjj, — > 77 and 7*7r° — > 7 in PJTOf, and Jain and 
Munczek's model passed this test very well. Other applications are also under investigation, 
and still many others are possible. However, for the full assessment of a model and for getting 
useful insights in how to improve it, it is also very interesting to see how it performs at the 
very edges of its applicability. Although Jain and Munczek's model is cleverly constructed 
so that it works well for most pseudoscalar and vector mesons below, above, and even on 
the mass scale of 77 and vj , the limitations of the ("improved" or "generalized" 0) ladder 
approximation employed by the model put the 77-77' system on such an "applicability edge" 
of this model - although not beyond it, contrary to what a pessimist could have concluded. 
This will be clarified below, where we analyze the rj-rj' system and its 77 decays in Jain 
and Munczek's model demonstrate the abilities and limitations of this model, and 

anticipate in which way it can be extended to improve further the description of 77 and 77' '. 

2. SOLVING THE CONSISTENTLY COUPLED SD AND BS EQUATIONS 

Dressed quark propagators Sf(q) for various flavors /, 
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Sf{q) = A f {q 2 )i-B f {q 2 ) , (f = u, d, s, ...) , (1) 

are obtained by solving the SD equation, which in the ladder approximation (i.e., with the 
true quark-gluon vertex replaced by the bare one, namely 7^A J /2) becomes 

Sj\p) = J- m f -ig*C F J -0^rS f (k)YG^(p - k) , (2) 

where mj is the bare mass of the quark flavor /, breaking the chiral symmetry explicitly, 
and Cf is the second Casimir invariant of the quark representation, here 4/3 for the case 
of the (halved) Gell-Mann matrices A J "/2 (j = 1,...,8) of SU(3) C . Neglecting ghosts, the 
product of the strong coupling constant g st and the Landau-gauge gluon propagator can be 



approximated by the Ansatz often described as the "Abelian approximation" |TT 



(k^k u \ 

As explained in the Introduction, the function G is given by the sum of the known 
perturbative part GW, an d the modeled nonperturbative part Gir. 

G{Q 2 ) = G UV {Q 2 ) + G IR (Q 2 ) , (Q 2 = -k 2 ) . (4) 

In Guy, we employ - following Ref. |J - the two-loop asymptotic expression for a st (Q 2 ): 

r tr* „ r ^ c " d r, , ,, '" |ln(x ° + ^ )1 1 

Guv ( Q)=iwCF —* ll + b ^ }. (5) 

where d = 12/(33 - 2N } ), b = 2(3 2 / (3\ = 2(19^/12 - 51/4)/(N f /3 - 11/2) 2 . As in §, we 
set the number of flavors Nf = 5, Aq C d = 228 MeV, and x = 10. We adopt the modeled 
Gir, together with its parameters a = (0.387 GeV)~ 4 and /i = (0.510 GeV) -2 , from Ref. ||: 

G IR (Q 2 ) =A-K 2 C F aQ 2 e-» Q \ (6) 

Solving (H) for the propagator functions Af(q 2 ) and Bf(q 2 ) also yields the constituent 
quark masses, defined (at q 2 = for definiteness) as A4 f = Bf(0)/Af(0) for the flavor /. 

The case mj = corresponds to the chiral limit, where the current quark mass rrif = 0, 
and where the constituent quark mass stems exclusively from dynamical chiral symmetry 
breaking (D^SB) 0. For u and d-quarks, the chiral limit is a very good approximation. 
Solving of © with m u = m d = leads to M u d = B u (0)/A u (0) = B d (0)/A d (0) = 356 MeV 
for the gluon propagator (|3|)-(|6|) with parameters quoted above and used in ||. 

When rHf ^ 0, the SD equation (0) must be regularized by a UV cutoff A an d rhf is 
in fact a cutoff-dependent quantity. We adopted the parameters of where (for A = 134 
GeV) the bare mass rrif (A 2 ) = 3.1 MeV - chosen to ultimately lead to the realistically 
massive pion - yields the light, non-strange isosymmetric constituent quark mass Ai u d = 375 
MeV, just 5% above its value in the chiral limit. For s-quarks, m/(A 2 ) is 73 MeV, giving 
us the strange quark constituent mass Ai s = B s (0)/A s (0) = 610 MeV |5|]. 
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In the chiral limit, solving of (|2[) with m/ = is already sufficient to give us the Goldstone 
pion bound-state vertex[] that is of zeroth order in the pion momentum p, 



A 3 ^ , 2 n , A 3 V2B f (q 



2\ 



TAq;p 2 = M^ = 0) = —T ff -(q;p' = 0) mf=0 = ^= l5 v ~ ' ^ ' ( 7 ) 

leading P^|,[7| to the famous result (Eqs. ( [ZED an d P8] ) below) for the 7r° — >■ 77 amplitude 
due to the Abelian Adler-Bell-Jackiw (ABJ), or axial, anomaly. 

Of course, for heavier qq composites one cannot circumvent solving the BS equation by 
invoking the chiral-limit (and the soft-limit, p^ — > 0) result (|7|). This is obvious when they 
contain c- or 6-quarks, for which the whole concept of the chiral limit is of course useless even 
qualitatively. When strange quarks are present, (|7p can be regarded only as an "exploratory" 
II expression, and it is useful for considering the chiral limit, since this limit is qualitatively 
meaningful for the s— quarks. Nonetheless, we need the quantitative predictions of Jain and 
Munczek's model for the ss pseudoscalar bound state, which is not physical, but enters as 
the heaviest component in the pseudoscalars 77 and rj', introduced in the next section. 

Therefore, we must obtain the bound-state vertex T sS by explicit solving of 

r„(?,J>) = iglC F J 0^l"S s (q' + l)T sS (q',P)Ss(q' ~ fh"^ - q') , (8) 

the homogeneous BS equation again in the ladder approximation, consistently with (fj). For 
pseudoscalar (P) quarkonia, the complete decomposition of the BS bound state vertex Y P 
in terms of the scalar functions Tf is: 

Tp(q,p) = 75 {r P (g,p) +]/rf (q,p) + ^(q,p) + ^J) T^q,p) } . (9) 

(The flavor structure is suppressed again. For neutral pseudoscalars, Tp is decomposed 



into //-components Ttj according to Eq. flTq) below.) The BS equation (jg) leads to a 



coupled set of integral equations for the functions Tf {i = 0,...,3), which we find to be 
most easily solved numerically in the Euclidean space by following the procedure of Jain 
and Munczek who formulate the problem in terms of the BS amplitudes XffiQiP) = 

S f (q + p/2)T fI (q,p)S f (q-p/2). 

In order to avoid the angular integration, we also adopt the momentum expansion (in the 
Chebyshev polynomials) of the four scalar functions appearing in the decomposition 

of the BS amplitudes. Ref. [f| often kept only the lowest order moment in the Chebyshev 
expansion, because they found it adequate for most of the meson spectrum. In contrast, 
while using the kernel and parameters of Ref. 0, we always retain all four functions in 



l ln (0), we explicitly included the (elsewhere suppressed) flavor factor A 3 /\/2, appropriate for ir°, 
to emphasize the change of our convention with respect to |3,10|: we now adopt the convention of 
Jain and Munczek's papers for the flavor factors, but not their conventional color factor of 
1/y/Nc. Hence, we have the additional factor of N c multiplying the integral in Eq. (2.8) of Ref. 
0], the formula which otherwise specifies our normalization. 
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solving of the BS equation (^, and the first two moments in the Chebyshev expansion. The 
accuracy of this procedure has recently received an independent confirmation - especially for 
presently interesting charge conjugation eigenstates - from Maris and Roberts [[11[]. In their 
study of 7r- and if-meson BS amplitudes, they employed both the Chebyshev decomposition, 
and straightforward multidimensional integration. Their comparison of these two techniques 
showed the very quick convergence of the Chebyshev expansion: in the case of equal quark 
and antiquark masses, such as in the pion, the zeroth and the first Chebyshev moment are 
enough for an accurate representation of the solution. Even for the kaon, still just one more 
is needed JTI] , in spite of the difference in the masses of its constituents. (Of course, the 
limitations of the ladder approximation would in the end lead to increasing difficulties if 
one of the fermion messes became much larger still, as recognized also by [||]. However, 
if the mass ratio of the constituents is not too large, various contributions beyond ladder 
approximation largely cancel out in the flavor-nonsinglet pseudoscalar, vector and axial 
channels fT^jnyTTfl , explaining the success of the ladder approximation in these channels.) 

Our procedure, already successfully used in P,|10|| for M Vc and M Vb , gives us M sS = 721 
MeV for the unphysical pseudoscalar ss bound state entering in the 7]-rj' system in the fashion 
discussed in the next section. Naturally, when we abandon the chiral limit approximation 
in (Q), we can also obtain the (isosymmetric) pion bound-state vertex T n o = T uU = 
replacing s — > u in (|8]). Although we stress that the chiral limit is an excellent approximation 
for many purposes in the case of pions, including the computation of tt° — > 77, it is also 
very important that the experimental 7r°-mass M„o = 135 MeV is reproduced || through 
(|) as M uU (= M d( j) with the small explicit chiral symmetry breaking, m u d(A 2 ) =3.1 MeV, 
corresponding to (isosymmetric) current u- and <i-quark masses m = 8.73 MeV, close to 
the empirical values extracted by current algebra. Such a small m cannot jeopardize the 
relevance of (0) for the computation of 7r° — > 77, as shown also by Ref. [T5[], which found 
(in an approach closely related to ours) that the amplitude decreased with respect to the 
analytic, chiral-limit axial anomaly result only by less than 1% when they introduced the 
non- vanishing but small u, d-quark mass m = 6.7 MeV. 



3. 77—77' COMPLEX AND ITS AXIAL-CURRENT DECAY CONSTANTS 

The SU(3)/ octet and singlet isospin zero states, rjs and r] , are in the gg-basis given by 

\ m ) = ^=(\uu) + \dd)-2\ss)) , (10) 

\Vo) = ^{\uu) + \dd) + \ss)) . (11) 

In our phenomenologically successful model choice ||, the flavor SU(3) /-symmetry is broken 
by the s-quark mass being realistically larger than the u, d-masses. Nevertheless, the isospin 
symmetry for u and <i-quarks is assumed exact throughout this paper. As is most commonly 
done, Eqs. ( |10D and (|Tl|) both employ the same quark basis states \ff) (/ = u,d,s) to 
define r] s and i]q. As pointed out by Gilman and Kauffman |16[ (following Chanowitz, their 



Ref. [8]), this usual procedure implicitly assumes nonet symmetry. However, it is ultimately 
broken by nonabelian ("gluon") axial anomaly, which will be discussed in Sec. |[ 
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r]s and t]q cannot be physical as they are not the mass eigenstates; that are their mixtures 
T) and rf\ 



\rj) = cos# \r] 8 ) - sin9 \rj ) , (12) 
\r]') = sin 6* \r]s) + cos 9 \r] ) . (13) 

The determination of the specific value that the mixing angle 9 should take, is a difficult 
issue which will be handled separately in Sec. |5|. We will keep our discussion general till we 
evaluate those of our results which are independent of the mixing and 9 - such as the decay 
constants of of the unmixed states r] s and r] - and point out those quantities for evaluation 
of which we need a concrete value of 9. 

For the light neutral pseudoscalar mesons P = tt°, i]g, 770, their axial-current decay con- 
stants fp = f n o, f m and f m , are defined by the matrix elements 

(0\m^Y^)\P(p)) = i^ p f P p^, (14) 

where if) = (u,d,s) is the fundamental representation of SU(3)/, while P = tt°, r] 8 ,r] si- 
multaneously has the meaning of the respective SU(3)/ indices 3,8,0. This picks out the 



diagonal (j = 3,8) SU(3)/ Gell-Mann matrices X j , and A = (^2/3)l 3 , in Eq. (0). 
The neutral pseudoscalars P are expressed through the quark basis states |//) by 




where the nonvanishing coefficients a p = (X p /y/2)ff for P = tt° are = — = \ = 
(X 3 /V2) u , whereas for ijg they are af = of = l/y/E = (A 8 /V2)n = (A 8 / v / 2) 22 , af = 
-2/y/E = (A 8 /V2) 3 3, and for V °, < = aP d ° = a? = (A /^)// = l/VS. 
The axial-current decay constants defined in flTJ]) can be expressed as 

fp= E ^T"///> (^ = ^ 3,77 8 ~8,and 77 ~0), (16) 

where we have for convenience introduced the auxiliary decay constant ffj, defined as the 
decay constant of the //-pseudoscalar bound state which has the mass Mfj and is described 
by the BS vertex Tjj(q,p), so that using the definitions of Bethe-Salpeter bound-state am- 
plitudes or vertices in the matrix elements (|i~4| ) as in, e.g., leads to 



/// = <^-rpf / (^itr [^S f {q + l)r ff -(q,p)S f ( q - \) 



(17) 



It turns out that this equation can also be applied for M^j = 0, as the limit exists. 

In the isospin limit, we get f n o = f uU = / d j = f ud = f w = 93.2 MeV in our chosen 
model ||. For the axial decay constant of the ss pseudoscalar bound state, we obtain 
f S s = 136.5MeV = lA7f T . (This factor with respect to f w is very reasonable and even 
expected, since the model also predicts the decay constant of the charged kaon f K + = 
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f uS = 114 MeV = 1.23/,,.) Eq. (pD then yields f m = 122.1 MeV and f m = 107.6 MeV. 
Note that f m = 1.31 f n , which is rather close to the result f m = 1.25 /V [[17]], obtained in 
the chiral perturbation theory (%PT). 

Evaluating the matrix elements of the pertinent mixtures yields the rj and rf decay 
constants 

cos 6> — sin 6> J f w + ( --^= cos# - sin#j / aS , (18) 



-y= cos 6> + sin 9 J f n + ( cos6> - sin6U / aS . (19) 



4. vr°, 778, 770 77 AND 77, 77' -> 77 PROCESSES 

The transition amplitudes for 77, 77' — > 77 can be obtained from the 77-transition amplitudes 
for t] 8 and t)q by forming the appropriate mixtures, in line with fllOP-fllBT). The 1] 8 ,1]Q ~ > 77 
amplitudes are in turn calculated in the same way as 7r°, t] c , r/ b — > 77 in . 

This means that we assume that these decays proceed through the triangle graph (de- 
picted in Fig. 1), and that we calculate the pertinent amplitudes 



T£*(Jfc, k') = e a ^ v k a k^T P {k 2 , k' 2 ) , (20) 
and the corresponding on-shell (k 2 = and k' 2 = 0) decay widths 

Vy(P^77) = ^M||T P (0,0)| 2 , (P = n°,r ] ,r ] ',...), (21) 
using the framework advocated by (for example) |l^,|^|T^J^,|l^] in the context of electromag- 



netic interactions of BS bound states, and often called the generalized impulse approximation 
(GIA) - e.g., by fi~5|j8[|. To evaluate the triangle graph, we therefore use the dressed quark 
propagator Sf(q), Eq. (^), and the pseudoscalar BS bound-state vertex Tp(q,p) instead 
of the bare 75 vertex. Another ingredient, crucial for GIA's ability to reproduce the cor- 
rect Abelian anomaly result, is employing an appropriately dressed electromagnetic vertex 
r^(q',q), which satisfies the vector Ward-Takahashi identity (WTI), 

(</ - q)^ f (q', q) = Sftf) - Sj l (q) , (/ = % d, s, ...). (22) 

Namely, assuming that photons couple to quarks through the bare vertex 7 M would be incon- 
sistent with our quark propagator, which, dynamically dressed through Eq. (|j), contains the 
momentum-dependent functions Af(q 2 ) and Bf(q 2 ). The bare vertex 7^ obviously violates 
(p2f), implying the nonconservation of the electromagnetic current and of the electric charge. 
Since solving the pertinent SD equation for the dressed quark-photon vertex is a difficult 
problem that has only recently begun to be addressed p0| , it is customary to use realistic 



Ansdtze. Following, e.g., |T5],|8],|7|,|iTfl , we choose the Ball-Chiu |2TJ vertex: 
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rjftf.g) = A{( q -, q ^ + ^^{^(g' 2 , 9 2 )^±£ - I^(g' 2 ,g 2 )} , (23) 

where H±(q' 2 ,q 2 ) = [Hf(q' 2 ) ± #/(g 2 )], for H = A or B. This Ansatz: - i) satisfies the 
WTI fl2~2|), - reduces to the bare vertex in the free- field limit as must be in perturbation 
theory, - Hi) has the same transformation properties under Lorentz transformations and 
charge conjugation as the bare vertex, - iv) has no kinematic singularities, and - v) does not 
introduce any new parameters as it is completely determined by the quark propagator (|I|). 
For the meson P whose flavor content is given by Eq. ([II]), GIA yields the amplitude 

Tj?(k, k') = ^ £ a p f Q) N c (-) J (0itr{r?(g -§,* + <?- \)S f {k + q-\) 



f=u,d,s 



x V f (k + q-Zq+ P -)S f {q + P -)Y ff {q, p)S f (q - |)} + (* «-> k' , (jl <-> u). (24) 

The coefficients of various flavor components \ff) in P = 7r°, r] 8 ,r] , are given below Eq. 
(|j~5l). Qf denotes the charge of the quark flavor /. The dependence on the flavor / has been 
indicated on the BS vertices, dressed propagators and electromagnetic vertices in the loop 
integral for each quark flavor. It is convenient to separate out the and Q 2 dependence by 
denoting each integral (times (—N c ) ) in fl24|) the "reduced 77-amplitude" T^. The "reduced 



scalar amplitude" Tjj for the flavor / is then 

ffi(k, k') = e a ^k a k'pf ff (k 2 , k' 2 ) . (25) 

A. 77-amplitudes and 77-decay constants 

Regardless of what the chiral-limit solutions for the propagator (|T|) and the bound-state 
vertex (^) are in detail, T^j(0, 0) can be evaluated analytically in the chiral (and soft) limit 
|7|JT^], which is perfectly adequate for f — u,d, i.e., for a Goldstone P = n°. There, 



N 

7>(0, 0) = T uu (0, 0) = T dd (0, 0) = , (26) 

to which we stick throughout. In terms of 

T P {k 2 ,k' 2 )^Y.^Q}T fI {k 2 ,k' 2 ), (27) 
/ 

this leads to the standard form of the successful axial-anomaly-resultQ for 7r° — > 77: 



We can also get - in the fashion of Ref. [19| — the anomalous "box" amplitude for 7* — > 7T7T7r. 
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T *° M =dk^ = 4^r (28) 

Note that this reproduction of the chiral limit relation between the 7r° — > 77 decay amplitude 
and the pion axial-current decay constant f n , is not dependent on the pion's internal struc- 
ture (or the interaction kernel that produces it) in any way J7|,|12| - and this is an important 
advantage of the coupled SD-BS approach over most other bound-state approaches, since the 
axial anomaly is on fundamental grounds known to be independent of the structure. (Those 
calculations of 7r° — > 77 which rely on the details of the hadronic structure - be it in the con- 
text of the BS equation without D^SB, nonrelativistic quarks, or otherwise - have problems 
to describe this decay accurately even when the model parameters are fine-tuned for that 
purpose; e.g., see p2|-p5| and references therein. The most successful of these model fits, 
Ref. |[24| , numerically obtains the width of 7.6 eV at the expense of fine-tuning constituent 
quark masses to unusually small values.) Of course, /„- itself is structure dependent. It is 
a calculated quantity in the SD-BS approach. Our model choice || successfully reproduces 
the experimental value of and this is obviously of utmost importance for the theoretical 
description of anomalous processes. 

The implications thereof for the r] 8 , r] and their mixtures r\ and rf are now clear, because 
those parts of their 77-decay amplitudes which stem from their uu and dd components are 
(just as in it ) accurately given by the Abelian anomaly (i.e., Eq. fl2"6])) for any interaction 
kernel which leads to the correct f n - be it the present one, or some improved one. 

In other words, Eq. (^) implies that the uncertainty (in the 77-amplitudes) due to 
modeling of the interaction kernel and the resulting bound state, is to a large extent cornered 
into the ss sector, since only T S g(0, 0) - the 77-decay amplitude of the relatively heavy ss- 
pseudoscalar - has to be evaluated numerically. From (f2~j|)-(p5|), we find numerically that 
in the model of §, f sS (0,0) = 0.62 f ufi (0, 0). 

The 7T° — > 77 decay amplitude 7X0,0) at any pion mass can be used as a definition 
of pionic 77- decay constant f_ n through 7X0,0) = l/^ 2 f n = {N c /2^2tt 2 : f n ) E/ af Q) . 
Eq. fl2"8p then reveals that f n = f n in the chiral limit, which result is well-known from 
the axial anomaly analysis. Although the chiral limit formula (p8|) can be applied without 
reservations only to pions, it is for historical reasons customary to write the amplitudes for 
V8,Vo ~^ 77 i n t ne same form as (p8|), defining thereby the 77-decay constants f m and f Vo : 



^ 0) ^?^^|^ (29) 



7X0, 0) = T ay Q) = h- ^ T «f>°\ (30) 



As pointed out by ||26|| , f m and f Vo are not a priori simply connected with the usual axial- 
current decay constants f m and f vo , in contradistinction to the pion case, where f n = f n 
because the chiral limit is such a good approximation for pions. 

Eqs. (|2?D-(|30"D reveal that in the present approach f V8 and f Vo are naturally expressed 
through f n (i.e., through T uU (0,0) = T d j(0, 0) evaluated in the chiral limit), and T s5 (0,0), 
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the 77-decay amplitude of the unphysical pseudoscalar ss bound state, calculated for non- 
vanishing m s . Our predictions for f vs and f m are thus: 

3 fn 

fm = 5-47r*y/2f„f aS (0,0)/N c ' (31) 



6 f 

fvo = -=-—^ . (32) 

Jm 5 + 27rV2AT ss -(0,0)/iV c 



Derivation of Eqs. (pi]) and (|32|) shows that irrespective of any specific model choice, any 
qq bound-state approach (such as our coupled SD-BS approach in conjunction with GIA) 
which has the merit of reproducing the anomalous 7r° — > 77 amplitude in the chiral limit, 
Eq. ( [ZED or (El)' should give the relations ( |3"T| ) and (^) for f m and f Vo , when pions are 
approximated by the chiral limit. The concrete numerical values of f m and f vo depend on 
what are the model predictions for f n and T S g. 

Since in the coupled SD-BS approach we can numerically evaluate T sS (0, 0) for arbitrary 



values of the s-quark mass, Eqs. (plf ) and (|32T ) give our predictions for the effects of the 
SU (3) / breaking on 77-decays in the rj-rj' system. In the SU (3) / limit (where T sS = T ua ) and 
the chiral limit applied also to s-quarks, we obviously recover f m = f n , but also f m = f n , 



since nonet symmetry in the sense of JTBJ is a starting assumption of ours. 

For our present model choice ||, where T sS (0, 0) = 0.62T uiZ (0, 0), Eqs. (|3l|) and (|32| ) give 

f m = 73.64 MeV = 0.797/^ , (33) 
f m = 98.58 MeV = 1.067/,, . (34) 



While this f vo agrees with the results of [P^7j^7| , there is a difference concerning / 



since f m > f n in %PT. This is important because the value of f n / f V8 has impact on the 
possible values of the 77-7/ mixing angle. We therefore devote the following subsection to 
the discussion of this result and the meaning of this difference. 

B. In SD-BS approach, f m < generally 



The /^-value (|33| ) is a result of a specific model. However, for the s-quark mass realistically 
heavier than the u, <i-quark masses, f V8 < f n holds in the coupled SD-BS approach generally, 
i.e., independently of chosen model details. To see this, let us start by noting that f Va < f n 
is equivalent to T V8 (0, 0) > T ff o(0, 0)/a/3, and since we can re- write Eq. ( p^) for r) 8 as 

^(0,0) = ^ + ^[^( 0j 0)-f ss -(0,0)] , (35) 

the inequality f V8 < f v is in our approach simply the consequence of the fact that the 
("reduced") 77-amplitude of the ss-pseudoscalar bound state, T sS , is smaller than the cor- 
responding non-strange 77-amplitude T rf( j (= T U u = T n o in the isosymmetric limit), for any 
realistic relationship between the non-strange and much larger strange quark masses. 

Only in the chiral limit (and close to it), subtle cancellations between the bound-state 
vertices, WTI-preserving qqj vertices and dynamically dressed propagators lead to the large 
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anomalous amplitude (|26|), or its slight modification (the size of which is controlled by 
Veltman- Sutherland theorem) for small u and d masses. Significantly away from the chiral 
limit, what happens is basically simple suppression of 7//(0, 0) by the large quark mass in 
the propagators in the triangle loop of Fig. 1. Essentials and generality of the suppression 
mechanism can be understood in basic terms in two (related) ways, through the simple free 
quark loop (QL) model and the Goldberger-Treiman (GT) relation. 

i) In a QL model (e.g., see and references therein), the strength of the Yukawa 
point couplings of the free quarks of the flavor / to the pseudoscalar P is given by the 
constant gj, and quarks have constant constituent masses hAf (in contradistinction to the 
momentum-dependent mass functions Aif(q 2 ) = Bf(q 2 )/Af(q 2 ) in our framework). Up 
to some arcsine-type dependence unessential here, each flavor / then contributes simply 
(gpff/Mf) Q} = (g / / M. /) Q 2 to the triangle-loop 77-amplitude p8[ . In the case of the 



strictly SU(3) j-symmetric coupling, the Yukawa couplings would be the same for all flavors, 
gf = g. The broken SU(3) /-symmetry implies that gf can differ for various flavors / - but 
not by much, so that relative strengths of the factors gf/A4f for various flavors is essentially 
determined by 1/M.f. Actually, this is what we find in our SD-BS framework, where the 
pseudoscalar bound-state vertices T ff are analogous to the coupling gf in the QL model, and 
gj/M.f is analogous to our "reduced" amplitude 7//(0, 0). Obviously, our approach allows 
for the flavor dependence of our BS Pqq vertices T^j, but due to the fact that the broken 
SU(3)/ is still an approximate symmetry, their variation with the breaking, given in terms 
of strange-to-nonstrange constituent mass ratio, is rather weak and cannot influence much 
the suppression occurring as the constituent mass in the denominator grows significantly. 
Hence, essentially the same mechanism is at work as in the QL model. That this parallel 
works very well, can be seen from the fact that the inverse of the strange-to-nonstrange 
constituent mass ratio in our SD-BS model, namely 1.63, quite accurately reproduces the 
suppression of the ss decay amplitude T sg (0,0) = 0.62T U „(0, 0), found numerically from 

ii) A related way to see the same effect is to apply the quark-level GT relation, gf/M.f = 
1////, to the pseudoscalars with the // quark content in the QL model. Then, roughly the 
same suppression factor occurs again, due to f sS = 1.47 f n . This is only roughly, since ss 
is further away from the chiral limit than uu and dd constituting the pion. Nevertheless, 
invoking the GT relation is in fact a very robust way to show that T sS (0, 0) < T uU (0, 0) must 
surely hold, even though s-quarks are much lighter than c- or 6-quarks (where the suppression 
is by orders of magnitude and D^SB is for s-quarks of importance similar to that for 
u, d-quarks. Precisely because the chiral limit makes sense for s-quarks qualitatively (as the 
pseudoscalars containing s-quarks, can still be considered pseudo-Goldstone bosons), the 
GT relation must continue to hold approximately in the ss-sector, regardless of any specific 
interaction kernel and of the resulting hadronic structure. T sg (0, 0) < T ua (0, 0) is therefore 
obligatory simply due to f sS > f n . 

The GT relation is useful also for demonstrating the robustness even of our model- 
dependent result on T S g(0, 0) - namely, that in spite of the model-dependence in the s-quark 
sector, the model kernel of our choice [§] should not lead to 77-amplitude T s5 (0, 0) exces- 
sively different from the ones which would result from an improved kernel. The usage of 
the GT relation at the quark level is especially transparent in the context of the simple 
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free quark loop model in which the GT relation, g//Mf = l////> ^ s necessary for repro- 
ducing the 77 anomaly amplitudes (|26|) and (|28D . In the context of the coupled SD-BS 
approach, with its dynamically dressed quarks and BS vertices, the GT relation for quarks 
and Goldstone bosons is given by the chiral-limit relation (^) (see, e.g., f29j). The chiral-limit 
analytic derivation of Eq. fl26|) from (|7|), with its subtle interplay and cancellations between 
the bound-state vertex, WTI-preserving qqj vertices and dynamically dressed propagators, 
transparently demonstrates the way the GT relation works for 77-decays in this context. 
For massive pions, the 77-amplitude must be evaluated numerically, but in fact changes 
very little, implying that the GT relation continue to hold very accurately. As just argued 
above, in ii), the GT relation should still hold as a rough approximation for the ss pseu- 
doscalar bound state. This, together with Eq. (|26|), implies that T sS ~ N c /(2\/2ir 2 f sS ). For 
f sS = 1.47/* which we obtained in the model |5[], this gives the GT relation-based estimate 
T sS ~ 0.68T U „. This is indeed in expected rough agreement with the accurate, numerically 
obtained prediction of the model §, that f sg (0,0) = 0.62T nS (0,0) = A^ c /(2v / 2tt 2 1.61/*). 
This shows that our result is quite reasonable. We remark that any model which is success- 
ful enough to reproduce empirical values of /* and fx+, should give a value for f sS close 
to ours, since anything very different from the estimate f sS ~ /* + 2(f K + — /*) would be 
unreasonable. Bound state descriptions that would be obtained by using kernels supposedly 
better than ours (improved beyond the ladder approximation by, say, including fully the 
gluon anomaly), must retain the good feature of agreeing approximately with the GT re- 
lation. This means that improving interaction kernels and, consequently, // bound states, 
would not change very much the 77-amplitudes with respect to our Tfj even in the s-quark 
sector. 



Regardless of any specific model realization of the coupled SD-BS approach, Eqs. (|5ID 
and Q52D with T sS (0,0) < T 7r o(0,0) imply the following bounds on f m and f vo . The equality 
holds when the chiral limit is applied to all three flavors, implying that /* is the upper bound 
for f V8 and lower bound for f Vo . As the s-quark mass grows, T S s(0,0) gradually diminishes, 
so that the lower bound for f V8 is 0.6/*, and the upper bound for f Vo is 1.2/*. 

Let us now address the meaning of the apparent contradiction between the results of the 
coupled SD-BS approach on f m j /* and the corresponding results of x?T, as well as possible 
ways to - at least in principle - overcome it. To this end, let us recall Pham's paper |3(| on 
gg-loop corrections S to the %PT result on f Va / f n . In the notation of our Eq. (p9[), he gets 

^(0,0) = A. (l-5)^M, (36) 

where the axial-current decay constant f Vs appears. From the standpoint of our approach 
and its Eq. ( |2"5| ) relating the 77-amplitudes of rj 8 and the chiral pion, / % /(l — 5) obviously 
corresponds to our f V8 . Can 6 resolve the discrepancy between our f m and f m = f m = 1.25/* 
of ;\TT fpTrl ? Pham's rough estimate is —0.28 < 8 < —0.19. This would practically reduce 



the %PT result from 1.25/* down to /*. In addition, his (following Ref. |31| ) result 



-!d + ln^r#) (37) 



3 V M 2 J \ A 2 

can be even more negative than —0.28, because there are large uncertainties in the value 
of his cutoff A (which can be lower than the lowest value of 1.2 GeV used by p0[), and in 
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his effective s-quark mass M. Besides, the current s-quark mass m s can be higher [32| than 
m s = 175 MeV used by Pham. Each of these possibilities would make 5 more negative. 
On the other hand, the ;\PT results |I7[], f m = f m on the equality of the axial and 77 
decay constants, and f V8 = (1 + ag)^ = 1.25 /V on SU(3)/ breaking effects in f m , are 
obtained in the one-loop approximation. Higher loops can introduce significant changes, 
including polynomial terms in meson masses, which may be both large with respect to 
chiral logarithms, and also not unambiguous [5^,17]. It is thus possible that as < 0.25 after 
all, which could make even easier for the correction factor (1 — 5) to reduce f vs , as defined 
by us in (F2D|), even below /„-. 

Large uncertainties in the quantities entering in the estimate for qq correction (|37"D, 
obviously leave plenty of room for usage of qq bound state models such as ours, which 
properly embed the chiral behavior of the underlying theory. Including meson loops in qq 
bound state approaches is a difficult task and implies going beyond the ladder approximation, 
but it would help further diminish the gap that Pham |JU| started closing from the side of 
%PT by incorporating into it the corrections due to quark loops. 

C. j], j]' —* 77 decay widths 



The 77, rj — > 77 amplitudes are given in terms of the 77-amplitudes of rjs and 770 as 

r„(0, 0) = cos 6T V8 (0,0) -sin TV (0,0) , (38) 
2V (0,0) = sin TV (0,0) + cos 9T m (0,0) . (39) 

Expressing TV (0,0) and TV (0,0) through the 77-decay constants f m ( EHf ) and f Vo (Effi), we 
arrive at the standard (e.g., see p6fl ) formulas for the rj and 7/ decay widths: 



W(r] 



77) 
77) 



a 1 M 6 

"em lvl r) 

64tt 3 3/2 

647T3 3/2 



/tt n 

-=r- cosy 
J 1 



fn . n 

sin t) 

J 1 



V8^ sin0 

JVo 

V8^ cos0 



1 2 



(40) 
(41) 



The version of (HO) and (Hip in which the axial-current decay constants f m and f Vo appear in 
place of f m and f Vo , requires a derivation where PCAC and soft meson technique is applied 
to r] and 7/ [p4| . These assumptions are impeccable for pions (leading to f n = f n ), but not 
for the 77-77' complex. In fact, the latter is quite dubious for the heavy 7/ ||34|| . We do not need 
and do not use these assumptions since we directly calculate the 773,770 — > 77 amplitudes, 
i.e., f m and f m . We also calculate f m and f vo independently of the 77-processes. 

For the values of f m and f m obtained in our model choice ||, namely (|33D and (^4|) , the 



best achievable consistency with the present overall fit |JIJ to the experimental widths, 
W exp (r] -> 77) = (0.46 ± 0.04)keV and W exp (7]' -» 77) = (4.26 ± 0.19)keV, (42) 



then occurs for = exp = —12.0° (obtained through x 2 minimization). Then 

W(r) -> 77) = 0.54 keV , 
W(rf -> 77) = 5.0 keV . 



(43) 
(44) 
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However, the present approach is capable of predicting the mixing angle 6, and it remains to 
be seen if the predicted 9 can be close to the angle favored by the experimental 77-widths. 

The issue of predicting 9 will be addressed in the next section. There is another mixing- 
independent quantity related to f m and f Vo , which we can predict before predicting 9. It is 
the i?-ratio, which is in fact measurable because it is the combination of 7r°, 7/ and rj widths: 



R 



W(t] — > 77) W(r]' — > 77) 



M 3 
77 



+ 



M 3 



Ml 



W(ir° 77) 



1 / f 2 f 2 

— I H. _|_ 8 1L 
Q { f 2 ° f 2 



(45) 



which is presently not known with satisfactory precision; [[36|] quotes R exv = 2.5±0.5(stat) ± 
0.5(syst). A more precise value of the .R-ratio ( f4"5| ) should come with DA$NE's operation 
at its higher energy a/s = 0.15 GeV, as this will enable the measurement of 77 — > rj' |p6|| . 
Since R is independent of 9, it will most cleanly test our predictions (|31"D, (|32|). Precisely 



determined i? exp can also help finding out whether (a) f m < f w , as follows in the coupled 
SD-BS approach from T sl (0,0) < T uU {0, 0), or (b) f V8 = f V8 > f n as in %PT The 
present model gives us T sS (0, 0) = 0.62T ua (0, 0), so that R = 2.87 is obtained, which is well 
within the error bars of the present experimental average fl3ql . (Taking the chiral limit also 
for s-quarks would give R = 3, which is the upper bound for (flop in the present approach. 



R = 3 is still consistent with R exp [^6| within the present experimental accuracy.) 



We should note: i) T sS is a quantity which can be especially practically used in conjunc- 
tion with a more accurate R exp to narrow down the choice of models suitable for describing 
rjs and 770 (and ultimately r\ and rj'). ii) Precise experimental determinations of R can help 
to find out if there are other admixtures \X) [e.g., gluonium \gg), r/(1295), 7/(1440) (former 
l), rjc, ... ] to 7/ and 7/', on top of the mixture of 7/ 8 (TO) and 7/ (|TT|) . We can see both and 



izj if we use our predictions for f Vs and /, 



'Al- 



and (H), in (H), yielding 



i? 



25 



2 
9 



T s5 (0,0) 



(46) 



which shows that in our approach R depends only on one variable^ model-dependent quan- 
tity: T sS . This is because f m (^Tj) and f m (P2] ) follow from the fact that, for any inter- 
action kernel and resulting propagator and bound state solutions, 7^.0 (0,0) = T nS (0,0) = 
A r c /(2\/ / 27r 2 / 7r ) in the chiral limit, and that this result remains an excellent approximation 
for realistic m u and rrtd leading to empirical M w . Therefore, important variations in our 



predictions for f m and f Vo , and thus R (^6|) , can come only from T sS (0,0). The accuracy 
of T s5 (0, 0) depends on the quality of the bound-state solution, but regardless of concrete 
model choices and results, the general inequality T uU (0, 0) > T s5 (0, 0) > enables Eq. ( |4~6|) to 
provide the bounds 3 > R > 25/9 = 2.777... . Hence, if experiments establish R < 25/9 by 
a significant amount, this will most probably indicate that in r\ and rj' there are admixtures 
{e.g., glueballs) to r/ 8 ( |T0"D and r/ ([□]) which are "inert" with respect to the interactions with 
photons, because this can lower the bound R > 25/9 most efficaciously We will be able to 
address this in more detail after the discussion of the mixing, presented in the next section. 



3 As clarified above, getting a reasonable /„- first is obligatory for applications to 77-processes. 
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5. COPING WITH MIXING OF ETAS IN COUPLED SD-BS APPROACH 



The mixing angle 9 is often inferred from the empirical 77 decay widths of 77 and rj'. This 
is how we established - in Sec. [| that # exp = —12.0° is the empirically preferred mixing 
angle for the values of f m and f Vo obtained in our model, namely (|33|) and (|34]). On the 
other hand, the angle 9 is predicted by diagonalizing the 77-7/ mass matrix evaluated in the 



r] 8 -r] basis - such as the one predicted by our SD-BS approach and given in Eq. ( |47D below. 
Obviously, for a satisfactory model description of the 77-77' complex, the latter procedure 
should give the mixing angle close to the angle # exp required by the 77, rj' — > 77 widths. 

In the 77-77' complex, subtleties arise from the interplay of the mixing due to the SU(3)j 
breaking with the Ua(1) gluon axial ABJ anomaly, which couples to the flavor- singlet 770 
and removes the nonet symmetry. (For a simple introduction, see Sec. 12.8 of || and Sees. 
III-3, VII-4 and X-3 of |26|). Namely, in the coupled SD-BS approach, where the states 
with good SU(3)j quantum numbers are constructed from the // bound states (/ = u, d, s) 

obtained in Sec. |2|, the eta (mass) 2 matrix M 2 in the rj 8 -rj basis (|T0|)-(|TT|) is given by 



M 



~M 2 8 


Mlo~ 




_M 2 8 
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3 V^ss 

4{mi 



{M% + \Ml) 
M%) 



UMl + Ml) 



(47) 



if we neglect the gluon anomaly for the moment. In agreement with other cases when the 
gluon ABJ anomaly is not included, or is turned off, as in N c — > 00 limit (e.g., [|37|,f26|1 ) , 
the diagonalization of (f|7|) yields an 77 degenerate with the pion, M 2 = M 2 , and without 
the ss component, 77 = -j^(uu + dd), whereas 77' is a pure ss pseudoscalar, with M 2 , = M 2 S . 
This happens at 9 = —54.74° and is obviously analogous to the "ideal" or "extreme" mixing 
which is known to be a very good approximation for the mixing of the vector mesons u and 
0. We can thus note in passing that the present approach works well for the mixing of to 
and 0. Nevertheless, this scenario is obviously catastrophic for the 77-77' system (the Ua(1) 
problem), so that gluon anomaly must be incorporated into our SD-BS framework. Doing 
this on the fundamental level represents a formidable task in any case, a task no-one has 
accomplished yet. Moreover, an interaction kernel in the ladder approximation, such as the 
simple gluon-exchange one that we have in the present model, is inadequate for this task even 
in principle. Namely, by definition it does not contain even the simplest annihilation graph of 
a quark- ant iquark pseudoscalar into two gluons (and their recombination into another quark- 
antiquark pair) contributing to the processes such as the one in Fig. 2. The contribution 
of the gluon ABJ anomaly operator e^^F^gF* to the 77o-mass, M 00 , therefore cannot be 
captured through a ladder kernel even in the roughest approximation (leaving alone the 
issue of non-perturbative gluon configurations such as instantons). 

Therefore, some additional ingredients or assumptions must be introduced into the 
present model in order to cope with the 77-7/ system. Since going beyond ladder approxima- 
tion is not within the scope of the present work, the following scheme is the most sensible 
at this level: note that there is a standard way (see, e.g., P,|26|) to account for the anomaly 
effect by parametrizing it through the term added to the 770 mass, since only this sin- 
glet combination ( |TT| ) is coupled to the gluon anomaly, so that only its mass is affected by 
it. This corrects the Ua(1) problem arising in the mass matrix evaluated with the nonet 
SU(3)j-states r] 8 and r] . Let us do the same in our mass matrix (^7|): 
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M 2 00 ^ ^-M% + M 2 W ) + X v . (48) 

Of course, parametrizing the effect of the gluon anomaly is far from actually calculating 
it unambiguously. In particular, the quantities we calculated for the 770 under the assumption 
of nonet symmetry, f m and f Vo , are in fact also affected by the coupling of the gluon anomaly 
to 770. However, due to the large N c arguments, it makes sense to break nonet symmetry 
only on the level of the mass- shift parameter A^ while keeping our r] Q built of the same // 
bound-state vertices as r/ 8 , to which gluon anomaly does not couple. This is because the 
gluon anomaly is in the large N c limit suppressed |37|j26|| as 1/N C , so that in our f m and f m 
calculated within the nonet scheme, only the contributions of the order 0(1/N C ) are missed. 
Our scheme is therefore a controlled approximation on the level of large N c arguments. 

The results obtained below for the mixing-dependent 77, 7/ — > 77 widths also turn out to 
be reasonable, providing an a posteriori justification for our scheme. In the light of large N c 
arguments, such reasonable results are not accidental and can be expected beforehand. 

Let us also note that our assumptions are in fact shared by many other approaches, 
explicitly or implicitly. E.g., Gilman and Kauffman [16| employ in their analysis nonet 



symmetry or broken version thereof, pointing out that it is at least implicitly assumed by 
all who use the quark basis not differentiating between quark states belonging to the singlet 
from those belonging to the octet. Moreover, imposing the nonet breaking via introducing 
the additional parameter X v is basically the same way in which nonet symmetry is broken in 
the chiral perturbation theory (%PT). In %PT, one faces the problem of how to incorporate 
r] , shifted upwards in mass by the gluon anomaly, into the scheme that should involve 



Goldstone pseudoscalar mesons. Bijnens, Bramon and Cornet p7[ comment on the problems 
encountered when working with this ninth state, but stick to what they did earlier p8" |, 
namely including rjo in their notation) "in a simple nonet-symmetry context". Their 
Ref. conveniently parametrized the nonet of (pseudoscalar) Goldstone bosons in terms 



of nine fields entering in the lowest order Lagrangian consistent with current algebra and 
explicit breaking by the quark masses, but the effect of the breaking of U(1)a is included 
only via an extra mass term for 770. This is justified if one relies on large N c arguments, since 
rjo is indeed a Goldstone boson in the limit N c — > 00 |39|,p6|], and then 770-mass is introduced 
as an extra parameter on top of that, which basically corresponds to our scheme. 



Precisely in the light of x?T, our result fl34|) for f m appears very reasonable in spite 



of missing the contributions of 0(1/N C ), supporting our relying on the nonet symmetry 
scheme. Namely, it is in excellent agreement with the results of x?T, being right between 
f V0 « l.lf n quoted by and f Vo = (1.04 ± 0.04)/, of 0. 



Finally, the robustness of our 77-amplitudes to kernel variations, resulting from the good 
chiral features of SD-BS approach (as explained in Sec. £|.B), also supports our scheme. 

We therefore pursue the procedure of removing the Ua(1) problem by lumping the effects 
of the gluon anomaly into a single r/o-mass shift parameter as in (fj8|). Then, with our 
result M sS = 0.721 GeV, and with the experimental pion mass M n o = 0.135 GeV (which the 
present SD-BS approach readily reproduces when the strict chiral limit is relaxed ||), and 



with the choice A,, = 1.165 GeV 2 , we get 9 = -12.7°, which is very close to 9 cxp = -12.0° 
favored by the empirical 77 widths. Moreover, we then reproduce the experimental value 
of 77-mass, M v = M^ xp = 0.547 GeV. Admittedly, the rj' mass is then somewhat too high, 
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M v > = 1.18 GeV. Of course, we could in principle pick such a value of A,,, that this other mass 
M v / would be reproduced, or still another value (\ v = 0.677 GeV 2 ) to reproduce empirical 
M 2 + M 2 ,. Naturally, M n would then be spoiled, but this is not the main reason why the 
two latter possibilities are disfavored in the present approach. The main problem is that 
they yield so negative values of the mixing ang le (9 = -21.4° when M 2 + M 2 , is fitted, and 
9 = —22.8° when M 2 , is fitted to experiment), that they are incompatible with the present 
approach. In the present model ||, so negative mixing angles obviously yield unacceptable 
77 widths, since the empirical 77 widths favor the mixing angle 9 exp = —12.0°, and this 
speaks in favor of the first possibility, = 1.165 GeV 2 , leading to 9 = —12.7° ~ 9 cxp . Still, 
this 9 exp = —12.0° is the consequence of the particular model choice || which led to the 
value (|33|) and ([34]) for f m and f m , respectively. Although we explained in the previous 
section why T sS (and consequently f m and f Vo ) must be relatively stable to model kernel 
variations, it is desirable to have a criterion which is even less mo del- dependent. And indeed, 
we do have a reason why the coupled SD-BS approach in general prefers the first procedure 
leading to larger values of \ v , and, consequently, less negative values of 9. Namely, it turns 
out that since in the coupled SD-BS approach f m < f n for any realistic value of strange 
quark mass, the consistency with the experimental rj, v[ — > 77 widths is possible in this 
approach only for mixing angles less negative than roughly —15°. This is easily seen, for 
example, in Fig. 1. of Ball et al. [0, where the values of f V8 , J fn consistent with experiment 
are given as a function of the mixing angle 9. [It does not matter that they in fact plotted 
/»78(o) //t an d n °t /»7 8 (o)//t- Namely, they used Eqs. pO|)-(|4lD for comparison with the 
experimental 77-widths, just with f Vsm /f n instead of f Vs{0) /fn, so that the experimental 
constraints displayed in their Fig. 1 apply to whatever ratios are used in these expressions.] 
On the other hand, the more negative values 9 < — 20° give good 77, rj' — ► 77 widths in 
conjunction with the ratio f m / f n = 1.25 obtained by [117] in x?T. However, our approach 
belongs among constituent quark ones. In the next section we discuss why considerably less 
negative angles, 9 « —14° ± 2 are natural for constituent quark approaches in general. 

The procedure leading to 9 = —12.7° is also corroborated by the results of some different 
approaches - most notably, by the results of the instanton liquid model, where one can 
actually calculate the gluon anomaly mass shift instead of parametrizing it. As Shuryak 
f42fl pointed out, the instanton-induced interaction leads simultaneously to both light pion 
and heavy rf - i.e., the dynamics provided by instantons can take care of the effects of the 
gluon axial anomaly and provide the light pseudoscalars as the Goldstone bosons of D%SB. 
While the instanton-induced interaction may therefore be the main candidate which in the 
future one may try to include in the interaction kernel of the coupled SD-BS equations, the 
results of Alkofer et al. E3| in the framework of the instanton liquid model have already 
indicated that such an inclusion could easily lead to a calculated X v similar to its present 
parametrized value. Namely, Alkofer et al. |^3| find that due to instantons, the U(1)a- 
anomalous contribution (2Nf/ f)N/V must be added to the rj-rj' mass matrix. This term, 
corresponding to our A^, also has the value very close to our A,, = 1.165 GeV 2 ; it is equal 
to approximately 1.1 GeV 2 for their standard instanton density N/V = 1 fm -4 and their 
pseudoscalar decay constant / = 91 MeV. (Number of flavors Nf = 3.) This gives them 
9 f=s —11.5°, M v ^ 0.527 GeV and ~ 1.172 GeV, which is very similar to our results. 
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A. Values of the mixing- dependent quantities 



Once the mixing angle 9 has been fixed, the predictions for the axial 77 and rj' decay constants 
are found from Eqs. (0)-(0). 6 = -12.7° implies = 112.6 MeV and f v > = 117.1 MeV. 



This agrees almost perfectly with Scadron's |0J estimates f v ~ 1.22/V and f v ' ~ 1-28/^ 
obtained from the GT relations at the quark level for the strange-to-nonstrange constituent 
mass ratio M. s /M. u d ~ 1.5 (and for 9 advocated by Scadron [44|, which is, interestingly, 
the same as our favored 9 = —12.7°.) However, these values are somewhat higher than the 
experimental values f% xp = 94 ± 7 MeV g5j or 79 ± 9 MeV g§ and f™ p = 89 ± 5 MeV 
or 96 ± 8 MeV j|6| , deduced (under certain theoretical assumptions discussed below) 



by CELLO |45[ and TPC/27 |^6[ collaborations from the Q 2 -dependence of their measured 
v(v')l*l transition form factors T r/ ( r? /)(0, — Q 2 ) (in our notation), where k' 2 = —Q 2 7^ is the 

The same TPC/27 reference f46 



momentum-squared of the spacelike off-shell photon 7* 
quotes also another pair of experimental values, f( ex P 2 ) = 91 ± 6 MeV and f^f xp ' 2) = 78 ± 5 
MeV, which were obtained from the experimental decay amplitudes into two on-shell photons 
under the assumption that one can write T^/^O, 0) = l/47r 2 /^( r? /) by analogy with the axial 
anomaly result (|28"D for the pion. However, because of the large s-quark mass, as well as 
the masses of r\ and 7/ which are, respectively, 4 and 7 times larger than the pion mass, this 
procedure can yield only a rough qualitative estimate. 



On the other hand, our value of f v is much closer not only to Scadron's |14| estimates 
and to the value f v = 114 MeV of an approach || somewhat related to ours, but also to the 
model-independent result of x?T, that f v = 1.02/ 7r (/ A -// 7r ) 4//3 f47j. For the experimental 
ratio f K /f n = 1.22 ±0.01, this gives f v = (1.3 ±0.05) f n = 120 ±5 MeV g7[, for which both 
CELLO H| and TPC/2 7 @ results are too low. 

The experimental values f^ xp and f^f p were extracted from the CELLO |R| and TPC/27 
|46| data on the transition form factors T^/^O, — Q 2 ) assuming that the pole mass ^M) 
parametrizing their fit to the data, can be identified with 27rv / 2/ r?( y). Then, the pole fits 
to the data could join smoothly (as Q 2 — > 00) the perturbative QCD prediction [48] for 



T v ( v >)(0, —Q )., i.e., the pole fits would then agree not only with the QCD asymptotic form 
l/Q 2 , but also with its coefficient. However, note that the values of f^ xp and /ff p quoted 
above, are all close to m p / (27r\/2) = 86.4 MeV, indicating that a connection with the vector- 
meson dominance interpretation (that ~ m p) |45|j46[] may indeed exist at the inves- 
tigated range of Q 2 . On the other hand, since Gasser and Leutwyler's model-independent 
calculation p7| , Scadron's |44| GT estimates, Burden et al. ||, and the present approach, 
all agree that f^ 1 ) should be noticeably larger than the extraction of f^ xp and f^f p from 
the transition form factors T r/ ( r) /)(0, — Q 2 ) probably cannot be done accurately at the ranges 



of Q 2 investigated so far. That this is indeed so, is indicated by the experimental value |35 
f n o = 84.1 ± 2.8MeV (w m p /(2nV2) again), extracted by the same method. The central 
value is 10% below well established f^ xp = 92.4 ±0.3 MeV. This cannot be explained by the 
small isospin violation, indicating that f^ xp and f^f p could have been underestimated too. 

Our predictions for the 77 and 77' two-photon widths are also totally fixed now, being 
given by our f m and f vo used in ([|(]) and (flTf) , without any additional parameters to adjust. 



Our preferred angle 9 = —12.7° leads to the predictions (displayed also in Table I.) 



W(r] ^77) = 0.561 keV. 



(49) 
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W(rf -> 77) = 4.913 keV . (50) 

These predictions are at first sight not very successful since, according to Table I, our best 
predictions overshoot the present |35|] experimental averages fl42|) for 77, rj' — > 77 by some 
20%. However, we should not be dissatisfied with these results because of the following: 
a) Ball et al. [^0| and, in effect, Review of Particle Properties itself (referring to the 



note on p. 1451 of |3l|), suggest that only the more recent data on 77,7/ — > 77 should be 



retained, whereby the presently "official" values (|42| ) are modified to |32|j40] 

W^iwin -> 77) = (0-510 ± 0.026) keV , (51) 
W e N x £ w ( V ' - 77) = (4.53 ± 0.59) keV , (52) 

and these experimental values agree much better with our predictions. 
b) We did not vary any model parameters, but used the parameters obtained from Ref.'s |5j 
broad fit to the meson spectrum and pseudoscalar decay constants. This fit did not include 
77—77' system in any way, so that everything we calculated for it are pure predictions. 
B. A side issue: speculations about other admixtures 

In the present approach, rjo and r/g (and consequently 77 and 7/) are constructed exclusively 
of the ground state pseudoscalar qq bound states. Nevertheless, it has been often speculated 
about additional admixtures, notably glueballs. Farrar jl9] points out that experiments 
appear to indicate that there is a glueball-like pseudoscalar which is much lighter than esti- 
mated by quenched lattice calculations, thus motivating us to speculate on the consequences 
of such admixtures. We do not have at this point the ambition to include such additional 
admixtures in our approach. However, we can anyway look into some of the consequences 
that such admixtures would have, by simply assuming that they were present besides the 
quarkonium r] and r] 8 as constructed in this paper. 

Take for example the simplest and most usual assumption |p6 , that only the SU(3)/ - 
singlet (|TI]) can be significantly modified in this way: 

COS iP — 

|r/ ) -> — t=-(|um) + \dd) + \ss)) + sin^|X) , (53) 
V3 

where <p is the new mixing angle, a new parameter expressing the assumed strength of the 
unspecified admixture \X) into rjo- 

If \X) is a state that does not couple to photons directly (e.g., gluonium \gg)), the results 
for 77-decays will be modified in a particularly simple way: in formulas fl30|) , (|40"D, ([41]) and 
(H), one should just replace l//o by cos </?/ fa. This can reduce R (|46]) strongly, as the 
largest term in Eq. (fBf), 25/9, would then be modified to 25/27 + cos 2 ip 50/27. 

We should also note that such an admixture (|53|) would help to fit the masses of both 
r) and r]' to their experimental values precisely - thanks to the new free parameter ip, of 
course. Eq. (^) modifies elements of the mass matrix to 

M| ^cos^(M 2 -il4), (54) 
M 2 - cos 2 V \(\m% + Ml) + X v , (55) 
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where = + sin 2 p> M\ takes the place that A,, has for p> = 0. If \X) is not a single state, 
but a mixture of various states, its mass Mx has a meaning of an effective mass. 

The experimental masses M v = 547 MeV and = 958 MeV, as well as the 77—77' mixing 
angle 9 = -17.1°, are then obtained for <p = 42.43° and X v = (0.873 GeV) 2 . Nevertheless, it 
turns out that the fit to the data is still not improved as much as one would expect when an 
additional free parameter is introduced, so that we did not detect indications for the need for 
an admixture of such states to what we have in the present model. For example, our i?-ratio 
then drops to R = 1.80. This is much further from the present central experimental value 
than R predicted by our approach without glueballs, but just in case that data from future 
precision measurements may strongly violate our bound on the i?-ratio, it is important to 
point out that - at least from the standpoint of our approach - such a violation would be 
a strong indication of the presence of some "inert" admixture, like gluonium. At present, 
however, the data are consistent with the bound R > 25/9 following generally from the SD- 
BS approach without gluonium admixture, and even favor the value R = 2.87 following from 
the present concrete model choice U without glueballs, over the value with the admixture 
quoted above. Moreover, the rj — > 77 width with the gluonium admixture improves only 
marginally, by 4%, while the 77' — > 77 width gets spoiled by more than a factor of 2. 

We therefore conclude that we found no indication that admixtures of glueballs, or 
other states with similar effects on 77-decays, would be favored by the present experimental 
data. Consequently, there is no strong motivation for enlarging the present framework by 
finding solutions for pseudoscalar glueballs and treating them on the same footing as our 
pseudoscalar qq bound states. (It is amusing that ip = 42.43° in conjunction with the 
vanishing gluon anomaly contribution, A^ = 0, implies Mx = 1.294 GeV - practically the 
same as the mass of 7/(1295). However, this can only be viewed as accidental at this point.) 
6. SUMMARY, DISCUSSION AND CONCLUSIONS 



The relativistically covariant constituent qq bound-state model H used here is consistent 
with current algebra because it incorporates the correct chiral symmetry behavior thanks 
to D%SB obtained in an, essentially, Nambu-Jona-Lasinio (NJL) fashion, but the model 
interaction is less schematic. Notably, when care is taken to preserve WTI of QED, it 
reproduces (in the chiral limit even analytically and independently of the internal meson 
structure) the Abelian axial anomaly results, which are otherwise notoriously difficult to 



reproduce in bound-state approaches (as illustrated by, e.g., |22[ and especially references 
therein). Observables such as meson masses, /„-, fx, fri, frf and 77-decay amplitudes can be 
calculated without additional parameters after an Ansatz has been made for the gluon prop- 
agator entering in the SD-BS equations, which are consistently coupled in the generalized (or 
improved) rainbow- ladder approximation (in the terminology of, e.g., or [jj). However, 
to avoid U(l) ^-problem in the 77-77' complex, we have to introduce an additional parameter, 
A^, representing the contribution of the gluon axial anomaly to the mass of 7/0, in analogy 
with the similar 7/0-mass parameter in the xPT Lagrangian in Ref. ||38|| , for example. Since 
the gluon anomaly contribution vanishes in the large N c limit as 1/N C , our qq bound-state 
pseudoscalar mesons behave in the N c — > 00 and chiral limits in the same way as those 
in xPT (e.g., see ||7]] or PE|]): as the strict chiral limit is approached for all three flavors, 
the SU(3)/ octet pseudoscalars including rj become massless Goldstone bosons, whereas the 
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//-mass is of order 1/N C since it is purely due to the gluon anomaly. In the N c — > oo limit 
with nonvanishing quark masses, the "ideal" mixing takes place so that rj consists of u, d 
quarks only and becomes degenerate with tt, whereas rf is the pure ss pseudoscalar. 

In our bound-state approach, f w , f m , f Vo , as well as f m , f vo and f v and f v i, are all 
calculated quantities, while most other theoretical frameworks treat at least one of them, 
f Vo , as a free parameter (fixed together with 9 from the experimental widths of rj, rf —>■ 77). 

Our prediction f m / f w = 1.31 agrees rather well with f m / f n = 1.25 of xPT [O. Nev- 



ertheless, this one-loop xPT calculation also lead to the identification of their axial-current 
and 77-decay constants, f m = f V8 , which differs from our results on r)$. More precisely, the 
observation that for realistic s-quark masses, T S s(0,0) < 7^0(0,0) always holds in the cou- 
pled SD-BS approach, leads to | f n < f m < f n and f n < f Vo < § fn- These inequalities hold 
irrespective of the model parameters or the quality of the interaction kernel. f m = f n = f m 
is realized in the chiral limit, whereas the opposite bounds are approached when the s-quark 
mass grows huge, leading to the decrease of T sl (0, 0)/T 7r o(0, 0). There is no disagreement with 
xPT regarding f vo , either concerning the general bound f n < f vo < § fn of the coupled SD- 
BS approach, or our result obtained using the concrete model of |J, namely f Vo / f n = 1.067. 
This agrees well with the values found in xPT []T7| , p7| . The apparent contradiction between 
the results of the coupled SD-BS approach on f w / f m , and the corresponding results of xPT 
was discussed in detail in Sec. [|.B. Let us now address the intimately related issue of 
different preferred mixing angles in these respective approaches. 

In conjunction with the updated experimental widths (|51| ) -(|52| ) , f m < f n implies that the 
coupled SD-BS approach is compatible with the mixing angles which are less negative than 
6 ps —15°. For our concrete model choice and the resulting values (p3|)-(p4[) of f m and 
f vo , the favored value of 9 is between the values accepted till mid-eighties, namely 6 ps —10° 
determined from the SU(3)/ breaking given by the Gell-Mann-Okubo mass formula, and 
the lowest of the values 9 G [-17°, -20°] favored nowadays j40|,pj|,p7|,p . 



In order to see that the mixing angles considerably less negative than those in %PT 
{9 ~ —20°) are a natural and expected prediction in a constituent approach such as ours, 
it is instructive to recall the paper of Bramon and Scadron j|l] where the mixing angle of 
9 = —14° ±2° follows from a rather exhaustive set of data if the SU(3) f breaking is taken into 
account in terms of the constituent quark mass ratio M. s / ' M. u d ~ 1-4 — 1.5. SU(3) ^-breaking 
ratios somewhere around this interval are considered realistic because they lead to good 
descriptions of many hadronic properties in numerous dynamical models; notably, close to 
this interval is also the ratio (~ 1.63) of the constituent masses Bf(0)/Af(0) generated by 
D^SB in Jain and Munczek's approach. Bramon and Scadron [f|l[] extracted their average 
9 = —14° ± 2° from the strong interaction tensor T — > PP decays, and the vector V — > 7P 
and pseudoscalar P — > 77 radiative decays. (And when extracted just from rj, rf — > 77 
pertinent here, and other SU (3) /-breaking-ratio-dependent radiative decays, their angle is 
even lower, —11° ± 2.4°.) They point out that more negative values 9 ~ —20° in the %PT 
framework are due to the way of implementing the SU(3) ^-breaking (through the values of 
the decay constants f V8 and f vo ), differing from that in the constituent-quark approaches. 

Now, our SU (3) /-breaking is fixed by Jain and Munczek's || choice of parameters, so 
that our calculated value of 9 varies only if we vary A,, which parametrizes the effects of the 
gluon anomaly. In the light of Bramon and Scadron's |41] observations discussed above, and 
the fact that that our SU(3) ^-breaking leads to the ratio of strange-to-nonstrange constituent 
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masses of 1.63, it is understandable and expected that our constituent approach should give 
reasonably good description of 77,?/ — > 77 for angles less negative than in ;^PT; i.e., it is no 
longer surprising that our preferred angle turned out to be 9 = —12.7°. However, it is not 
only that these values are the preferred ones in our presently chosen model || because they 
are more empirically successful than other values. In addition to that, since in the coupled 
SD-BS approach f Vs < f n rather generally (for any realistic value of strange quark mass), 
the consistency with the experimental rj, rj' — > 77 widths is possible - in our approach - only 
for mixing angles less negative than roughly —15°, as already pointed out above. 

That all this is in qualitative agreement with what was known from relatively simple- 



minded constituent-quark models even before the analysis of ||41[|, can be seen, e.g., from 



Zielihski's review [[50] on radiative decays of mesons. He observed that in the scenarios that 
related apparent suppressions of radiative decays of strange mesons to a larger mass of the s- 
quark, a significant suppression of the annihilation amplitude of ss pairs into two photons was 
also expected, and with the latter suppression of order 0.5 relative to annihilation amplitudes 



of non-strange quarks (relevant, e.g., to the model of ]5T]]), the two-photon widths of both 



rj and rf could (however roughly) best be described with 9 ~ —11°. Remembering the 
limitations on mutually consistent 9 and f V8(0) , we see that our values of 9 and fit 

with our third element, T sl (0, 0) = 0.62T u{i (0, 0), into a logical scheme which is consistent 
with the behavior of the approaches similar to ours. Zielihski [[)tj also discussed how 9 was 
much more negative (~ —20°) in chiral theories, but pointed out that the determination of 
the pseudoscalar mixing angle was model dependent, and a clean-cut choice among various 
schemes was rather difficult to establish. Our discussion, and the results of, e.g., Bramon 
and Scadron [41|, Pham [30] and Ball et al. [40], shows that this assessment still holds, but 
also that there has been some progress in narrowing the interval of possible mixing angles. In 
this connection, recall the observation of |4(|], that newer experimental input [our Eqs. (|51~D 
and (0)] reduces the mixing angle even more than Pham [[30] realized, to 9 = —(17 ± 2)°. 
This is not any more so far away from our preferred 9 (especially considering that the value 
of the correction 5, Eq. (|3"TD, can be even more negative than Pham's values |50[). If various 
approaches succeed in including physical mechanisms they have been missing so far, their 
predictions for 9 will probably tend to a unique value. Likewise will be with f m and f vo . 
In view of |4l] , |3l3| , f4~0"|| , this final value at which 9 will settle, may well be roughly in between 
the values favored nowadays by and by quark model approaches such as ours. Thus, 
9 ~ —14° to —17° maybe encompasses the final result. An f m which would be rather close to 



the chiral limit value f m = f n , because the chiral-loop contributions would be - like in |30 
- to some extent (over) canceled by some other contributions (like our bound-state strange 
mass-breaking effects), would agree better with such a 9 ~ —14° to —17°. In our approach, 
the physical mechanisms which are now absent, are those corresponding to loops in %PT. 
Including them obviously implies substantial enlargements beyond the present framework. 
However, this also holds for others - e.g., in x?T one might pose the question what would 
the effects of higher loops and vector mesons be. At present, no approach can claim to have 
all the relevant physics included, and therefore the ultimate values for 9 and f m . 

The present experimental value of the i?-ratio (f|^) is described reasonably well by our 
approach. What if more precise measurements (e.g., at DA$NE |36[]) constrain R eX p below 
25/9? A strong violation of this bound (say, R exp < 2.5) would indicate that important 
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admixtures other than r)s and rjo are present in 77 and rj'. If the violation is not that strong, 
the following possibility is also viable: some of the values in the interval 2.5 < R exp < 25/9 
can be satisfied by f V8 and f m predicted by %PT. Hence, such smaller violation of our bound 
can also mean that the prediction of xPT, that f n < f m , is favored over our prediction. 
This would indicate that in the case of the rj-rj' complex, the ladder-approximated SD-BS 
approach makes larger error by neglecting meson loops than, e.g., in the case of the charge 
pion form factor calculated in in the context of SD equations, where the contribution of 
meson loops was estimated to be much smaller than that of the quark core ||52|| . 

The quantities dependent on the 77-77' mixing, namely axial-current decay constants f v and 
f v >, masses and 77-decay widths of 77-7/, are satisfactorily close to data (or other theoretical 
predictions such as %PT) considering that - except for parametrizing the mass shift due to 
coupling of t]q to nonabelian axial anomaly - we did not do any parameter fitting, but used 
the parameters obtained from Jain and Munczek's || broad fit to the meson spectrum and 
decay constants. We conclude that their model [||] again performed well. 

Since the coupled SD-BS approach is - due to the key role of D^SB - akin to the NJL 
model conceptually, the progress we made is best illustrated through the comparison with 
the analysis of the tt°,i] — > 77 decays and properties of the pion, kaon and 77, performed in 
a NJL model (extended to include three flavors and the 't Hooft determinantal instanton- 



induced interaction) in ]53] and in parts of (7/ was not treated in |53 



For the choice of model parameters preferred by Takizawa, Oka, and Nemoto |53|-|55 



the experimental amplitude for 77 — > 77 is reproduced, but the 77-mass is 7% below the 
experimental value. The mixing angle is 6 = —1.25°, showing that their ?7(1 ^-breaking 
is stronger than in our approach (not to mention the one in ^PT), forcing their 77 to be 
an almost pure Their kaon decay constant fx = 96.6 MeV is 15% below the observed 
one. Accordingly, the predicted 77 decay constant, f v ~ is uncomfortably far from what 
the model-independent result of xPT g7|, = 1.02f n (f K / f n )^ 3 , gives when the empirical 
JkI U is plugged in. 

While our results compare rather favorably with the above, the best examples of ad- 
vantages both in the conceptual consistency and in the quantitative details which coupled 
SD-BS approach has with respect to the NJL-model, are P — > 77 decays. Namely, we must 
criticize the seemingly successful reproduction of the anomalous 7r°, 77 — > 77 amplitudes by 
|53^p5[. In contradistinction to the coupled SD-BS approach with nonlocal interactions - 
in particular Jain and Munczek's model, where the UV cutoff is either not needed (in the 
chiral limit or practically infinite compared to the relevant hadronic scales - the NJL 
approach contains a low cutoff. In spite of this, Refs. ||5T3l-j5"5]j leave the convergent integrals 
unregulated, because the triangle diagram reproduces the anomalous tt° — > 77 amplitude 
(p8|) only if there is no NJL cutoff |53|j56| -p8| - else, an underestimate of, typically, 20% 
occurs ||57|| . While Refs. |^3j|55l claim the improvement of the rj — > 77 amplitude and width 
with respect to the earlier treatment of Bernard et al. [|59| , the consistent viewpoint is that 
of [^J : once the cutoff is introduced, the effective theory is defined and should not be altered 
for the purpose of calculating various quantities. In such effective theories, the missing part 
of the anomalous amplitude - lost due to the cutoff - should be found in additional diagrams 
|56| which contribute since the cutoff cannot be let to infinity. I.e., in the class of models 
employing only local interactions and therefore needing a low cutoff, a simple incorpora- 
tion of the anomaly is not possible p5]-p^], in contrast to the coupled SD-BS approaches 
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employing also nonlocal interactions and thus not having such a cutoff. 

Ref. || is another approach to qq substructure incorporating D^SB, and it is even closer 
to us than the NJL model. The interaction used in || is nonlocal, like ours, allowing the 
generation of momentum-dependent dynamical mass and BS-vertices, so that there are no 
problems with a low cutoff like in the NJL model. The mixing angle they favor, 6 ~ +5, 
results from its treatment as an external parameter on which the mass and other properties 
of rj depend. However, their axial current decay constant f v = 114 MeV is close to ours. 

Extending the treatment of the 7r°7* — ► 7 transition form factor of ||[10| to the 77(7/) 7* — > 
7 transition form factors is presently under investigation pU . 
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TABLES 



p 


W(P -► 77) 


W exp (P 77 : 






7.7 


7.74±0.56 


not applicable 


V 


0.56 x 10+ 3 


(0.46 ± 0.04) x 10+ 3 


(0.51 ± 0.026) x 10+ 3 


rf 


4.9 x 10+ 3 


(4.26 ±0.19) x 10+ 3 


(4.53 ± 0.59) x 10+ 3 



TABLE I. Comparison of the calculated 77 decay widths (in eV) of vr°,ry and r( with their 
average experimental widths, as well as the experimental widths W^ew obtained when only more 
recent measurements are taken into account. The widths are calculated using the empirical masses 
in the phase-space factors in conjunction with calculated amplitudes. The tabulated r\ and rf 
calculated widths correspond to the case when their mixing adjusts the mass of r] to its empirical 
mass. 
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FIGURE CAPTIONS 



Fig. 1: The diagram for P — * 77 decays (P = 7r°, 77, 77', ...). Within the scheme of generalized 
impulse approximation, the propagators and vertices are dressed. 

Fig. 2: The annihilation of the // bound-state vertex T^j into two gluons and their recom- 
bination into the quark- ant iquark bound state consisting of possibly different flavors 

/'/'■ 
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